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' We derive an exact general formalism that expresses the eigenvector 
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and the eigenvalue dynamics as a set of coupled equations of motion in 
terms of the matrix elements dynamics. Combined with an appropri- 
ate model Hamiltonian, these equations are used to investigate the effect 



(N 

> . 

■r^ , of the presence of a discrete symmetry in the level curvature distribu- 

o 

tion. An explanation of the unexpected behavior of the data regarding 
^ . frequencies of acoustic vibrations of quartz block is provided. 

C ■ I. INTRODUCTION 

The usefulness of the study of statistical properties of eigenvalues and eigenvectors 

X 

' of quantum systems has already been demonstrated in many areas of physics. A 
lot can be learned, especially about symmetries, by just employing the appropriate 
statistics. It has also become clear that these statistics follow universal patterns that 
can be modelled by probability distributions extracted from an ensemble of random 
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Hamiltonians of the same class of the underlying symmetry of the system under study 
This has been a field of intense investigation over the last two decades 0. These 
activities have concentrated their effort on what we can call the "statics" of the problem, 
in which stationary Hamiltonians are considered. More recently, however, the interest 
has also been directed to the dynamical aspects of the same question. 

The "dynamics" consist in considering a given Hamiltonian as a function of a param- 
eter (representing "time"). The statistical properties that characterize the evolution 
are then studied as the parameter is varied. Only evolutions that preserve the sym- 
metry class of the Hamiltonian are considered. Several measures have been introduced 
to investigate this kind of evolution. One of the most used ones is the probability 
distribution of the level curvature, which can be thought of as "acceleration" as it is 
defined in terms of the second derivative with respect to the parameter. These dis- 
tributions measure correlations among the set of eigenvalues. Another measure that 
is commonly used is the two-point correlation function between first derivatives ("ve- 
locities"). Given some generic level, this two-point correlation function is obtained by 
calculating the "velocity" at two different values of the parameter P|. Measures have 
also been considered to characterize the evolution of the eigenvectors 0]. 

These studies started with Wilkinson's pioneering work that investigated the de- 
pendence of the eigenvalues of a fully chaotic billiard as a function of its shape 0. 
The plot of the trajectories of levels as a function of the parameter that controls the 
shape, exhibits a typical pattern of avoided crossings. A measure of these is provided 
by the curvature of the trajectory. There is now an analytical evidence that, in the 
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fully chaotic regime, the curvatures, after an appropriate rescaling, follow an univer- 
sal simple distribution. The tail of this distribution has been investigated, and an 
asymptotic dependence inversely proportional to the third power of the curvature was 
established, for fully chaotic systems that are time reversal invariant and thus governed 
by the Gaussian Orthogonal Ensemble (GOE) The expression 

P{k) 



2(1 + A;2)t 



was then proposed for the entire domain of the curvature k. Finally, it has been proved 
that this function gives the exact distribution of curvatures, in the case of random 
matrix ensembles 0. The power | in the denominator is the GOE value of ^^^^^ 
with /5 = 1,2 and 4 for the GOE, Gaussian Unitary Ensemble (GUE) and Gaussian 
Symplectic Ensemble (GSE), respectively. 

Recently, the difficult task of checking experimentally this prediction was under- 
taken by the experimental group at the Center for Chaos and Turbulence of the Niels 
Bohr Institute They studied the dependence on the temperature, the external 
parameter in the system, of the spectrum of frequencies of quartz blocks. In previous 
investigations , they have found that the spectra of frequencies of quartz blocks obey 
statistical models based on random matrix theories. The dynamics of the frequencies, 
as a function of the temperature, was therefore measured for a quartz block whose 
static statistical properties were already previously established. The data obtained, 
however, have shown a deviation from the above expected curvature distribution. This 
deviation, although slight, is significant and not yet completely understood. We are 



going to show, in the present paper, that the data can in fact be fully understood if 
one requires the average curvature to be equal to one, as is implied by the universal 
distribution, Eq. (|l]). 

So far, all studies of parametric correlations have been concentrated on the fully 
chaotic regime when the system statistics are well described by the Gaussian ensembles 
of Random Matrix Theory (RMT), in particular, the Orthogonal Ensemble (GOE), 
if there is time-reversal invariance. The partially chaotic situation has been little 
investigated. We intend here to provide the first systematic discussion of this situation. 
In section II, we develop the formalism and the model we employ, and in section III, 
we present the numerical results and discussion. We verify that, at the GOE limit, the 
above universal expression for the level curvature distribution is obtained. As some 
degree of symmetry is introduced, it is found that the distribution becomes narrower. 
However, as the symmetry is progressively introduced, the distributions return to the 
universal function, if the average of curvatures is imposed to be one. 

II. THE FORMALISM AND THE MODEL 

We shall now derive a set of equations that describe simultaneously the dynamics of 
the energy levels and of the eigenvector components of a Hamiltonian H . These equa- 
tions contain the equations of motion of the matrix elements of if, whose dependence 
on the parameter, t, representing the "time", is supposed to be given. Our starting 
point is the general matrix equation 
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H = UHdU\ (2) 

where H is an A^x real symmetric matrix, Ho is the diagonal matrix constructed with 
the N eigenvalues, and U is the unitary matrix whose columns are the N eigenvectors. 
Assuming that Ho and U also depend on the parameter t, differentiating Eq. (0) with 
respect to t we get 

H= U Hd U HdU^ + UHd U\ (3) 

where the derivative is indicated by a dot. Multiplying (0) by from the left, by U 
from the right and defining the matrix S = U'^ U= — U"^ U we obtain the equation of 
motion 

Hd=[Hd,S] + P (4) 

where the matrix P = W H U was introduced. On the other hand, we find for P the 
conjugate evolution equation 

p = [p, S] + U^HU (5) 

By choosing a particular model, i.e., the dependence of the matrix element on the 
parameter t, these equations can be employed in several contexts. They can be used, for 
example, to construct an alternative method of matrix diagonalization, or, by requiring 
the matrix elements to satisfy appropriate Langevin equations, they lead to Dyson's 



Brownian motion model [0. Here, we concentrate on the simple model given by 

H = Hicost + H2smt , (6) 



where Hi and H2 are a couple of fixed, i.e., parameter independent, random matrices 
taken from the same matrix ensemble, and t is the parameter. If in Hi and H2 are 
taken from a Gaussian ensemble, the evolution will preserve the probability distribu- 
tion, so that H will remain in the same ensemble. With this choice, becomes 

P=[P,S]-Hn (7) 

The pair of coupled equations, (^) and (|^, have the explicit solution 

Hd (t) = (t) [Hd (0) cos t + P (0) sin t] U (t) 

and 

P (t) = (t) [-Hd (0) sint + P (0) cost] U (t) 

where U {t) is the solution of the equation U= U S given by 

U{t) = Texp [ S{T)dT 
Jo 

with T being the time-ordering operator. 

To implement this solution numerically a basis has to be chosen to express the 
eigencomponents. Since our objective is to investigate the curvature, a quantity related 
to the behavior of the eigenvalues as the external parameter is varied, it is convenient 
to use the instantaneous Hamiltonian eigenstates as basis vectors. In this case, from 
the diagonal part of (^) we derive 

Ek = Pkk, (8) 
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and using in (H) the relation Ski = Pm/ {Ek — Ei) (obtained from the off-diagonal part 



of (ID) we can then derive the equations 

TV 
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m=l,m^kA 



(9) 



l,m^ky 

and 
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Pkk=-Ek+ y: e^z^- 

m=l,m7=fe 

This set of coupled equations is one of the main results of this paper. All calculations 
will be based on it. Thus the "accelerations", i.e. the levels' curvature, are just given 
by (0). 

Regarding the random matrix ensemble, we shall work with a Gaussian ensemble 
that interpolates between one GOE and two decoupled GOE's. This ensemble has 
been already employed with a very satisfactory result in the analysis of data relative to 
symmetry breaking |[TT|JT^ in nuclear [|l^] and acoustic systems 0]. It can be defined 



by the following operator equation [jl4| 



H = PH'^'^'^P + QH''''''Q + A {PH'^O'^Q + QH''''''P) , (li; 



M 



where P=^Pj, Q = l — P and Pj =| i >< i \, i = 1, . . . , N are projection operators. 



i=l 



< A < 1 is the parameter that controls the transition, and H'^'^^ denotes a GOE 
matrix whose elements follow a joint probability distribution given by 

P oc exp [-atr (i^^'^^)'] , (12) 

with a being an arbitrary scaling parameter. With the above definitions, A = 1 corre- 
sponds to the GOE case, while A = corresponds to block diagonal random matrices. 



made up of two GOE matrices of sizes M x M and (A^ — M) x {N — M) . 

Regarding this ensemble, it is important to stress that A is not the more convenient 
parameter to work with, since the transition is also dependent on the matrix size, N. 
Independence on the dimension is obtained by introducing the scaled parameter 

£ = ^/NX. (13) 



III. NUMERICAL RESULTS AND DISCUSSION 

Before presenting the results, we discuss the rescaling variables necessary to extract 
a universal behavior. First, we have to unfold the spectrum, that is, we work with a 
new spectrum generated by the transformation 

xi= dEp{E) ioT I = 1,...,N , (14) 



where p(-E') is the averaged level density. Without loss of generality, we consider in 
the calculation only the symmetric situation, = 2M, in which the matrices are 
decomposed into blocks of equal size. In this case, the average density is given by the 



Wigner's semicircle law |16]. With an appropriate scaling that guarantees the correct 



value of the second moment of the eigenvalue, the level density is given by 



In the Fig. 1, we show the nice fit obtained with this expression when compared with 
the numerical values of p (E) generated within the two coupled GOE's ensemble alluded 
to above. 
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Second, we need some normalization of the accelerations. This is a controversial 
issue that requires some discussion. On one hand, it has been proposed that the 



parameter t should be replaced by a new dimensionless parameter r related to t by |15 

|=VM. (16) 

where the average of the velocity is made over the whole set of eigenvalues or, equiv- 
alently, over the ensemble. The level curvature is then defined in terms of these new 
scaled variables as 

where x= p (E) E and x= p{E)E + ^ (^E^ ^ . 

On the other hand, the universal curvature distribution, Eq. (P, implies that 
(I A; I) = 1. It is not at all clear that the scaled curvatures given by equation ([T7|) will 
satisfy this condition. Thus, we imposed the normalization 

with K given by (|T7|) . Our calculations have shown that this last step is necessary in 
order to get stable results, i. e., independent of the subset of levels of the spectra over 
which the statistics is performed. 

The behavior of the distribution, Eq. (|I|), for large curvatures can be traced to the 
level spacing distribution. In fact, large curvatures can be considered, approximately, 
as inversely proportional to the small level spacing s. Thus if we assume s ccl/k and 
use the fact that, in the GOE case, P (s) is linear in s, we obtain 
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~ k-^ , (19) 

as predicted by (|l|). As a consequence, as symmetry is introduced by decreasing the 
parameter A, one would expect a reduction on the probabihty of large curvatures with 
the distribution becoming narrower. We shall see that this indeed happens. 

Our main results are presented in Fig. 2, where the curvature distributions were 
calculated for four values of the scaled parameter (|13D. The figures show that, at the two 
extreme situations, namely, in the one GOE limit, and in the two fully decoupled GOE's 
limit, the curvatures distribute themselves according to the universal distribution. We 
would expect this kind of behavior in the latter limit, since the levels in each block 
become completely independent of the levels in the other. Thus, their trajectories can 
cross freely. As the discussion in the preceding paragraph predicted, the distributions 
are narrower in the intermediate region. We stress that these results are strongly 
dependent on the renormalization ([T^). 

Turning now to the question of the behavior (distributions wider than the universal) 
presented by the data of Ref. P], one possible explanation would be that the curvatures 
do not average to one. To check this point we have fitted the data with the distribution 

P (K) = ^ (20) 

27 [1 + (K/^f] ^ 

in which the average curvature, 7 =< \K\ >, is treated as a free parameter. The best 
fit, obtained with 7 = 1.27 ±0.01, is displayed in the figure 3. This excellent fit makes 
this explanation plausible. 

In conclusion, we have investigated the effect of the symmetry breaking on the level 
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Pik) ~ P(s) 



ds 



dk 



curvature distribution using a random matrix ensemble that allows for the transition 
from one GOE to two decoupled GOE's. We have also provided an explanation for 
the discrepancy of the data regarding the temperature dependence of frequencies of 
acoustic vibrations of quartz blocks. 
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Figure Captions 

FIG. 1. Densitiy of levels: comparison of the calculated density (histogram) with the 
semi-circle law ( |T5|) (solid line). The calculation corresponds to matrices of dimension 
N = 100, and £ = 0.32. 

FIG. 2. Level curvature distributions: comparison of the calculated histograms with the 
theoretical prediction (|l]) (solid line). The calculations correspond to matrices of dimension 
N = 100, and for the values of e indicated in the figure. 

FIG. 3. Fitting of the data (crosses) of Ref. Q with the parametrized distribution (pO|) 
(solid line). The best fit (x^ = 0.00004) was obtained with 7 = 1.27±0.01. The dotted curve 
corresponds to the universal distribution. 
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